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Abstract 

The shock-induced combustion of methane-air mix- 
tures in hypersonic flows is investigated using a new re- 
action mechanism consisting of 19 reacting species and 
52 elementary reactions. This reduced model is derived 
from a full kinetic mechanism via the Detailed Reduc- 
tion technique. Zero-dimensional computations of several 
shock-tube experiments are presented first. The reaction 
mechanism is then combined with a fully implicit Navier- 
Stokes CFD code to conduct numerical simulations of 
two-dimensional and axisymmetric shock-induced com- 
bustion experiments of stoichiometric methane-air mix- 
tures at a Mach number of M = 6.61. Applications to 
the ram accelerator concept are also presented. 

Introduction 

One of the main obstacles in the computation of high- 
speed flows with hydrocarbon combustion has been the 
lack of reliable reaction mechanisms that are reasonably 
simple and yet still reproduce experimental observations 
over a wide range of conditions. Detailed kinetic mecha- 
nisms have been developed for the simplest fuel molecules. 
Oxidation models for methane, for example, typically 
consist of 100-250 elementary reactions, and more sig- 
nificantly, they include from 28 to 50 or more species^ 
Since the computational cost associated with a given re- 
action mechanism depends primarily on the number of 
species included, the use of these models would be pro- 
hibitively expensive when combined with existing two-or 
three-dimensional CFD codes. 

As a result, previous computations of high-speed 
methane-air flows, aimed at studying detonation waves 
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and the ram accelerator concept^, have used either 
simplified quasi-global models’^-®, or mechanisms com- 
posed of elementary reactions whose rate parameters are 
collected from literature recommendations®. Although 
quasi-global models have the potential to reproduce ac- 
curately some flame properties such as flammability lim- 
its, combustion temperature and burned gas composition, 
they generally cannot describe accurately the chemical 
structure of the flame itself. In shock-induced combustion 
problems this will translate in an incorrect prediction of 
the induction zone. In addition, the parameter ranges 
over which these simplified mechanisms can accurately 
be applied are rather narrow. A second disadvantage of 
quasi-global models, such as those developed by West- 
brook and Dryer^, is their sensitivity to the numerical 
models used. These models, which combine a single re- 
action of fuel and oxidizer with a detailed mechanism for 
a CO - jS ’2 - O 2 system, were developed mainly based 
on flame speed data. As a result, they depend not only 
on the rate parameters, but also on thermodynamic and 
transport properties which may be treated somewhat dif- 
ferently in other codes. Therefore, for use in other codes, 
it is generally necessary to calibrate the pre-exponential 
factor in each individual reaction. 

On the other hand, mechanisms that are simply a col- 
lection of elementary reactions suffer from the following 
principle: 

‘‘A mechanism composed of reactions with best avail- 
able rate parameters individually is incapable of quanti- 
tative predictions when taken as a whole^”. 

For these mechanisms, in which the individual reactions 
may have been validated under conditions of chemical 
isolation, when combined together they usually do not 
reproduce experimental observations accurately. 

There is therefore a need for a reasonably simple and re- 
liable kinetic mechanism for methane-air combustion. In 
the present paper we introduce one such model which con- 
sists of 19 reacting species and 52 elementary reactions. 
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Although this model is still quite complex and requires 
significant computational power, it has the advantage of 
giving accurate predictions over a wide range of flow con- 
ditions. Its use in practical hypervelocity applications is 
demonstrated in this paper. The combustion model is 
combined with a fully implicit CFD code to numerically 
simulate expansion-tube experiments conducted at the 
French-German ISL Institute^® aimed at studying ram 
accelerator related combustion phenomena. Results are 
presented for superdetonative flows of methaine-air mix- 
tures around various bodies, and are compared with ex- 
perimental data and with the computations performed 
by Soetrisno, Imlay and Roberts®. Applications to an ax- 
isymmetric ram accelerator projectile are also presented. 

Numerical Formulation 


Rabinowitz^. Detailed Reduction is a systematic method 
of reducing large reaction networks while maintaining the 
accurate prediction of selected combustion characteris- 
tics. Usually a kinetic mechanism is developed to predict 
a range of combustion characteristics (ignition, species 
profiles, flame velocities, pollutant emission, etc.) many 
of which are not important for a particular application. 
For detonation calculations the fate of species that do 
not contribute significantly to heat release or ignition is 
unimportant. To test the contribution of various species 
and reactions two criteria were developed: 

\Kj\<eK\Kref\ (3) 

and 

< egQma. (4) 


Governing Equations 

The computations are conducted using the Navier- 
Stokes equations for two-dimensional or axisymmetric 
flow, in which the global continuity equation is replaced 
by all the species continuity equations. They can be ex- 
pressed in the following conservation form for a gas con- 
taining n species and in general curvilinear coordinates 


£Q_^a(F — d{G G„) ^ -yy (1) 


where 
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The equations describe two-dimensional flow if j = 0 
and axisymmetric flow if j = 1. The variables are the 
density of the ith species pi, with p = Pi’> the ve- 
locity components u and v, and the total energy per unit 
volume e. F and G are the inviscid flux vectors in the 
^ and rj directions respectively. Similarly, F„ and G^ 
are the viscous fluxes. The terms H and are the ax- 
isymmetric source terms, and W is the chemical source 
term. A detailed description of all the terms appearing 
in Eq. (1) can be found in Yungster^^. 


Reaction Model 


where, Kj is the rate of reaction Kref is the rate of a 
reference reaction (usually the rate limiting reaction, in 
this case the reaction of H O 2 ^ OH -f O), AHj is 
the enthalpy of reaction j, Qmax is the maximum heat 
release per unit time of any reaction, and €k and eq are 
parameters much smaller than unity. The first inequality 
tests the contribution of each reaction to chain branching 
(and hence ignition delay) and the second to heat release. 
Reactions whose rates, both forward and reverse, satisfy 
equations 3 and 4 are removed from the mechanism. 

A series of zero dimensional calculations were per- 
formed with values of ck and eg around 0.1 for 11 
shock-tube test cases of premixed methane oxidation. It 
was found that the complete 33 species, 149 reaction 
mechanism^ could be reduced to 19 species and 52 reac- 
tions while still maintaining good accuracy. The results 
are shown in Table 1. The computed flame velocity for a 
stoichiometric, 1 atm methane flame is 39.2 cm/s, in close 
agreement with both experimental value (40 cm/s) and 
the value obtained with the full mechanism (39.5 cm/s). 
The reactions for the reduced mechanism and their rate 
coeflicients are listed in Table 2 in the appendix. 

The procedure for calculating the chemical source term 
TY is ais follows: 

The chemical equation for a general elementary reac- 
tion j in a gas mixture containing n species can be written 
as 


( 5 ) 

3=1 3=1 

where u' • and v'' ■ are the stoichiometric coefficients of 
reaction j. The rate of change of the concentration of 
species i in reaction j, denoted as i/ij^ is given by 


The methane oxidation mechanism used in this study ^ „ 

was developed via the technique of Detailed Reduction^^ vij = n y*'" - n (6) 
starting from the full mechanism of Frenklach, Wang and 
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Table 1: Ignition Delay (microseconds) 


Fc = aj€ + (1 - 0'j)e 


(13) 


IcT 

Experiment 

Full* 

Reduced’" 

1 

268 

236 (-12) 

245 (-9) 

2 

46 

49 (7) 

51 (11) 

3 

18 

22 (22) 

23 (28) 

4 

436 

497 (14) 

530 (22) 

5 

99 

101 (2) 

107 (8) 

6 

99 

97 (-2) 

99 (0) 

7 

1512 

1361 (-10) 

1378 (-9) 

8 

550 

538 ^2) 

546 (-1) 

9 

226 

249 (11) 

253 (12) 

10 

205 

187 (-8) 

192 (-6) 

11 

210 

241 (15) 

255 (21) 


tSee Ref. 1 for details of cases ICl through ICll. 


’"Percent deviation in parentheses. 


where yi is the concentration of species Yi {yi = pi/ Mi). 
The total rate of change of the concentration of the ith 
species caused by all of the chemical reactions is 

yi = Yly'd ( 7 ) 

j 

The components Wi of the chemical source term W are 
given by 


Wi = ViMi (8) 

The forward rate coeflBicients of pressure-independent 

reactions are calculated from 

= (9) 

Note that the dissociation reactions of and 

C 2 HQ (reactions 33, 48 and 49) are in the fall-off region 
and require special treatment for pressure-dependent rate 
coefficients. The pressure dependence of these reactions 
is treated based on the Troe-Golden formalism^ 



( 10 ) 

= AfT^fe-^riT 

( 11 ) 

" Kr 

( 12 ) 


where and K"p are the low and high pressure lim- 
iting rate coefficients and [M] denotes the concentration 
of the third body. Note that M has been written within 
parentheses in Table 2 to indicate that it should not be 
included in Eq. 6. 


l + (10gP.)2 

The forward rate coefficients of pressure-dependent re- 
actions are then given by 

Kf,^=KrY^/cXt (15) 

The rate coefficients of the reverse reactions were de- 
termined via equilibrium constants 


(16) 

where Kcj is the equilibrium constant for the jth reac- 
tion, and is given by 


Kc,i = 


exp{- Vsii^JUT) 


where 


(17) 


n 

^ V3) 1/5 

3=1 



for reactants 
for products 


(18) 


Here, R is the universal gas constant, and is the 
species standard state Gibbs free energy per mole, which 
is determined from fourth-order polynomials of temper- 
ature. (The other thermodynamic and transport prop- 
erties are also computed using polynomial functions of 
temperature^^). 


Numerical Method 


The system of equations (1) is solved using a fully im- 
plicit finite difference CFD code^^. It employs an itera- 
tive method that is based on the LU-SSOR implicit fac- 
torization scheme^^, and Yee’s second order total vari- 
ation diminishing (TVD) differencing scheme^^. In the 
present study, a symmetric TVD scheme with a minmod 
type limiter is used. The viscous terms are evaluated us- 
ing standard central differences. The full Jacobian of the 
chemical source term is used, leading to a preconditioner 
matrix of size n x n that has to be inverted at every grid 
point. The structure of the LU-SSOR algorithm allows 
full vectorization of the CFD code, including the matrix 
inversion process which is done using Grout’s algorithm 
(without pivoting). The scheme is suitable for steady- 
state calculations. Further details about the algorithm 
can be found in Yungster^^. 
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Results 

Shock-Induced Combustion 

Numerical simulations of two sets of expansion-tube 
experiments conducted by Srulijes, Smeets and Seiler^^ 
at the French-German ISL Institute are presented. Both 
cases considered a stoichiometric methane-air mixture at 
a static pressure, poo = 0.51 bar, static temperature. 
Too = 295°Ff, and a superdetonative velocity of Uoo = 
2330 m/s (Mach number M = 6.61). The Chapman- 
Jouguet detonation speed of the gas mixture under these 
conditions was computed^^ to be about D = 1800 m/s. 

The first set of experiments consisted of cylindrical steel 
rods placed perpendicularly to the flow inside a combus- 
tion test chanaber. For a fixed flow velocity, the diameter 
of the rods were progressively reduced up to a lower limit 
beyond which ignition did not occur any more. Two pres- 
sure transducers mounted on the tube wall, one placed 
upstream and the other downstream of the rod, recorded 
the pressure history. The results of this experiment in- 
dicated a sharp ignition onset between a rod diameter 
of 3 mm and 4 mm. For a rod diameter of 7 mm, the 
pressure trace clearly showed the presence of combustion. 



Fig. 1. Temperature contours for a rod diameter of 
d = 1 mm. 

Figures 1-3 show the computational results for three 
cylindrical rods having diameters of 1 mm, 3 mm, and 7 
mm respectively. A 91 x 91 grid was used, and the flow 
was assumed to be inviscid. The results show the advance 
of the combustion zone, initially constrained to a narrow 
region near the surface of the smallest rod. towards the 
bow shock as the diameter of the body is increased. The 


increase in the amount of combustion as a function of rod 
size appears to be a continuous process. Therefore, it is 
probable that the experimental setup of Ref 10 was un- 
able to detect the partial combustion predicted in Figs 1 
or 2 because the pressure rise was too small or occurred 
downstream of the location of the pressure transducer. 



Fig. 2. Temperature contours for a rod diameter of 
d = 3 mm. 



Fig. 3. Temperature contours for a rod diameter of 
d = 7 mm. 
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The same computations were conducted by Soetrisno 
et. al.® using the quasi-global combustion model of West- 
brook and Dryer®. Their computations predict essentially 
a fully coupled shock-deflagration wave for the same three 
rod diameters studied. These results show the inability of 
the global models in correctly predicting induction times. 
The computation of Soetrisno et. al.® for a rod diameter 
of 0.5 mm is comparable to the result presented in Fig 3 
of the present work for a 7 mm diameter cylinder. Ne- 
glecting the fact that a somewhat coarser grid was used 
in Ref 8 (72 x 65), this indicates that the global model is 
underpredicting the induction time by approximately an 
order of magnitude. 

Figures 4-5 show the pressure and temperature varia- 
tion along the stagnation streamline for the three cases. 
For comparison, the nonreacting solution is also plotted. 
The shock standoff distance increases with the amount of 
heat release. The pressure drop behind the shock wave 
is due to the heat release at subsonic speeds. Figure 6 
shows the mole fraction distribution along the stagnation 
streamline for the 3 mm cylinder. Note that behind the 
shock, thermal decomposition of CH^ produces signifi- 
cant amounts of larger hydrocarbons and other radicals 
within the induction zone. These species then quickly 
disappear during the heat release process. 

These simulations required between 2500 and 3500 it- 
erations (using a maximum CFL number of « 8) and 
between 7 and 10 hours of CPU time on a CRAY C90. 
The code achieved 175 MFLOPS during execution, of 
which virtually all were performed in the vector units. 
The program, however, appeared to have memory bank 
conflicts with this number of species. (The same code 
with an 9-species, 18-reaction model achieved 250 

MFLOPS). A more efficient method for referencing mem- 
ory is currently being implemented. 

The second experiment consisted of a blunt cylinder of 
d = 7 mm diameter, having its axis aligned with the flow. 
The mixture, Mach number, and free-stream conditions 
were identical to the previous case. Laminar, adiabatic 
flow alculations were performed on a two-block grid hav- 
ing 75 X 150 and 75 x 91 points. Temperature contours for 
this case are shown in Fig. 7. The shock and combustion 
front are separated by a small induction zone. This in- 
duction zone widens downstream due to the weakening of 
the bow shock caused by the expansion waves emanating 
from the cylinder shoulder. 

The same calculation was also conducted by Soetrisno 
et. al.® using several combustion models. Their results 
indicate also a decoupling of the shock and combustion 
fronts downstream of the shoulder, however, along the 
stagnation region, the shock and combustion fronts ap- 
pear to be fully coupled. 



Fig. 4. Pressure distribution along the stagnation 
streamline. 
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Fig. 5. Temperature distribution along the stagnation 
streamline. 
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Fig. 6. Concentration profiles along the stagnation 
streamline; d = 3 mm. 



Fig. 7. Temperature contours for a blunt cylinder. 
M=6.61, stoichiometric 

Figure 8 shows the pressure and temperature distri- 
bution along the stagnation streamline obtained in the 
present work, compared with the results of Soetrisno et. 
al. based on the quasi-global model and a detailed model 
consisting of 13 species and 19 elementary reactions. All 
computations give similar pressure distribution with the 
exception that the von-Neumann spike predicted in the 
present calculation is not observed in the calculations of 
Soetrisno et. aJ. The temperature distribution shows sig- 
nificant differences between the two combustion models 
used in Ref. 8, with the global model showing a better 
agreement with the present work. Also, the shock stand- 
off distance was slightly larger for the global model than 
that obtained in the present work. 

A quantitative comparison with the interferometry flow 
visualization of this case presented in Ref 6 could only 
be done with respect to the shock stand-off distance. A 
steady combustion process was reported, with a constant 
shock stand-off distance Xsh- The experimental value of 
Xah/d « 0.29 is in good agreement with the present re- 
sults and with the full model of Soetrisno et. al. 

Applications to the Ram Accelerator Concept 

The ram accelerator is a chemical propulsion method 
for accelerating projectiles to very high speeds®’^^. In this 
device, a shaped projectile is accelerated inside a tube 
filled with a premixed gaseous fuel/oxidizer mixture. In 
the high-speed modes of the ram accelerator, ignition is 
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achieved by means of a detonation wave or other forms of 
shock-induced combustion. Several methods for generat- 
ing the detonation wave have been proposed. In the sim- 
plest method, a series of shock waves reflected from the 
tube and the projectile surface increase the temperature 
of the mixture until the ignition temperature is reached 
at a designed location. The energy released will then es- 
tablish a detonation wave or a shock-deflagration wave, 
depending on the mixture composition, pressure and tube 
size. 




x/d 


Fig. 8 . Pressure and temperature variation along 
stagnation streamline. Cylinder diameter, d=7 
mm. Soetrisno^- detailed 13-species, 19-step model; 
Soetrisno^- quasi-global model. 


In a second method, the forward cone angle of the pro- 
jectile is kept small to reduce drag and prevent premature 
combustion, and a detonation wave is generated aero- 
dynamically by inserting a sudden, short axisymmetric 
ramp^^ with a relatively steep angle, or even a forward 
facing step as proposed by Rom and Avital^® in a sim- 
ilar External Propulsion Accelerator concept. Inviscid 
simulations of these oblique detonation drive concepts, 
using simple one step kinetic models, were conducted by 
BogdanofF and Brackett^*^ and Tivanov and Rom^®. Ac- 
curate prediction of the ignition delay is important for 
scaling studies, and viscous eflFects cannot in general be 
neglected since the boundary layer can interact with the 
detonation wave and with other shocks to strongly modify 
the inviscid picture. 

In this paper, we present flrst a laminar flow investiga- 
tion of the generation of a detonation wave by means of a 
double ramp using the newly developed methane oxida- 
tion mechanism. Then, we incorporate this concept in a 
ram accelerator configuration. 

Figure 9 shows a schematic of the wave structure for 
nonreacting and reacting flows. In the nonreacting flow 
case, the two shock waves 5 i and S 2 intersect at point 
0. From the intersection point O, there is a transmitted 
shock wave sz and a weaker wave which, at hypersonic 
speeds, is always an expansion wave. Regions 3 and 4 in 
Fig 9a are separated by a slip line sl\. The temperature 
in region 4 is always higher than that of region 3 (except 
for the boundary layer). This presents the interesting 
possibility that, for given free-stream conditions, ignition 
be achieved in region 4 but not in region 3. This has the 
advantage that the detonation wave does not interact di- 
rectly with the body. However, some combustion can still 
take place along the boundary layer. This situation, in 
which region 3 acts as a buffer zone between the detona- 
tion wave and the body surface, is illustrated in Fig 9b. 
Behind the transmitted shock S 3 , there is an induction 
zone 4, at the end of which energy release becomes sig- 
nificant and generates a set of deflagration waves through 
which there is a smooth rise in pressure and temperature. 
These deflagration waves converge into the shock S 3 , mak- 
ing it steeper until a new oblique detonation wave and a 
second slip line are formed. An additional compression 
wave S 4 is needed to equalize the pressures on the two 
sides of the slip line sl 2 - The flow structure in zones 4, 5 
and 6 is similar to that found in a basic detonation on a 
wedge^^. 

Figures 10 and 11 show temperature contours of a nu- 
merical simulation conducted on a double ramp, with 
0i = 18® and O 2 = 36®. The mixture considered is 
CH 4 -f 4 O 2 + 15.04iV2 (equivalence ratio <t> = 0.5). The 
free-stream pressure and temperature are poo = 1 atm. 
Too = 300® AT, and the Mach number is M = 7.5. The 
computations assumed laminar flow and a constant wall 
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temperature = 600®iT. The length of the flow domain 
is 40 cm. Besides the flow structure already described, 
note that combustion occurs in the boundary layer from 
the corner of the double ramp and downstream. As a 
result of the reacting boundary layer, a large separation 
bubble is established at the corner. 




Fig. 9. Schematic of the wave structure on a double 
ramp for reacting and nonreacting flow, s- shock wave; 
si- slip hne; E- Prandtl Meyer expansion wave; df- defla- 
gration waves; ODW- oblique detonation wave. 

Figure 12 shows the pressure and temperature varia- 
tion along the gridlines k = 32, k = 45, and k = 65. For 
gridline k = 32, the nonreacting plot shows the pressure 
and temperature jumps across the two shocks si and 52 , 
followed by the expansion across E. The reacting solu- 
tion along this same gridline shows in addition a jump 
caused by the separation shock and, behind the expan- 
sion, a pressure increase due to combustion. Note that 
the temperature remains near the nonreacting level ex- 


cept when this gridline crosses the boundary layer and 
a temperature rise is observed. Along the A; = 45 grid- 
line the pressure and temperature show the jump across 
sz followed by a short smooth rise due to the deflagra- 
tion waves and a final jump caused by the detonation 
wave. Note that the temperature decreases as the grid- 
line crosses the slip line sli. In the nonreacting flow case, 
the temperature in region 4 (Fig 9a) is T 4 « ISOO^iiT, 
high enough to ignite the mixture, while the temperature 
in region 3 is only T 3 « llOO^if, too low for ignition. 



Fig. 10. Temperature contours for nonreacting, M=7.5 
flow past a double ramp. 61 = 18®, ^2 = 36®. 



H 40 cm J 


Fig. 11. Temperature contours for reacting, M=7.5 flow 
past a double ramp. 0i = 18®, 62 = 36®. 
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Fig. 12. Pressure and temperature distribution along 
three ^ = const gridlines (indicated in Fig. 10). NR- 
nonreacting flow; R- reacting flow. Reference length L = 
5 cm. 


try is based on a 90 mm tube diameter similar to the 
ISL ram accelerator, and is shown in Fig. 13 . The pro- 
jectile considered has a 15® nose cone half angle and a 
30® axisymmetric ramp. The overall length of the pro- 
jectile is 18.95 cm. Although previous inviscid analy- 
ses have placed the ramp downstream of the projectile 
shoulder^^’^*^, viscous calculations indicated very strong 
shock/boundary layer interactions that strongly modified 
the ideal inviscid flowfield. The projectile configuration 
considered in the present work is similar to that proposed 
in Refs. 16 and 18. 

Computations were conducted on a 195 x 70 grid for 
both inviscid and laminar flow. Figures 14 (inviscid) 
and 15 (laminar) show temperature contours (top half) 
and CO 2 mass fraction contours (bottom half) for a 
CH 4 , -h 4O2 + I5.O47V2 mixture at poo = 1 atm, Too — 
300®iir, and a Mach number of M = 9. Under this con- 
ditions, combustion was initiated by the reflected shock 
from the tube wall. Steeper cone angles or a scale up 
of the ram accelerator is required in order for the trans- 
mitted shock to ignite the mixture in this case. The adi- 
abatic laminar result (Fig. 15) shows combustion along 
the boundary layer, and a large separation region behind 
the projectile shoulder caused by the interaction between 
the reflected detonation and the reacting boundary layer. 
Note also that the separated boundary layer affects the 
manner in which the expansion waves emanating from the 
projectile shoulder interact with the transmitted shock 
and with the reflected detonation. The pressure distribu- 
tion along the projectile surface and tube wall is shown in 
Fig. 16. The high pressure on the back of the projectile 
caused by the combustion process can be clearly iden- 
tified. In the laminar calculation, a small pressure rise 
occurs immediately behind the projectile shoulder due to 
the separated boundary layer. 

Conclusions 

A new reaction mechanism for methane-air combustion 
was presented. Although this model is still quite com- 
plex and requires significant computational power, it can 
be accurately apphed over a wide range of flow param- 
eters. Its use in practical hypersonic flow computations 
was demonstrated. For specific situations in which some 
species and/or reactions are presumed to be unimportant, 
the present model could serve as a benchmark against 
which simpler models could be compared. 


The plot for gridline fc = 65 shows the jump across 
the oblique detonation wave. Note that the pressure at 
the outflow boundary is basically uniform (as it should 
be) but that there are large variations in the temperature 
profile. 

The final set of calculations is conducted on an ax- 
isymmetric ram accelerator configuration. The geome- 
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accelerator barrel 


Fig. 13. Schematic of ram accelerator. Grid 195 x 70. 
Dimensions are in mm. 
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Fig. 14. Temperature contours (top half) and CO 2 
mass fraction contours (bottom half). Inviscid M = 9 
flow. 



Fig. 15. Temperature contours (top half) and CO 2 
mass fraction contours (bottom half). Laminar M = 9 
flow. 


Fig. 16. Pressure distribution along the projectile sur- 
face and tube wall; M=9. Reference length L = 2.3684 
cm. 
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Appendix 


Table 2: Methane-air reaction mechanismt 



Reaction 

A 

0 

b 

1 

H + 02^0H + 0 

1.59 X 10^'^ 

8491.28 

-0.927 

2 

0 + H2^0H + H 

3.87 X 10^ 

3151.16 

2.70 

3 

0H + H2^ H 2 O + H 

2.16 X 108 

1725.92 

1.51 

4 

OH + OH ^0 + H 2 O 

2.10 X 10» 

-199.65 

1.40 

5 

H + H + M^H2 + M 

6.40 X lOi’^ 

0.0 

-1.0 

6 

H + OH + M^ H 2 O + M 

8.40 X 10=^^ 

0.0 

-2.0 

7 

fl- + O 2 + M ^ HO 2 + M 

7.00 X lO^’^ 

0.0 

-0.8 

8 

HO 2 + H^OH + OH 

1.50 X 10^^ 

505.15 

0.0 

9 

HO 2 H ^ H 2 + O 2 

2.50 X 10^=* 

348.79 

0.0 

10 

HO 2 + 0 ^ O 2 + OF 

2.00 X 10^2 

0.0 

0.0 

11 

HO 2 + OF ^ F 2 O + O 2 

6.02 X 10^^ 

0.0 

0.0 

12 

H 2 O 2 "f A/ ^ OH -f- OH -}- Af 

1.00 X 10^^ 

22851.89 

0.0 

13 

CO + OH^ CO 2 + F 

1.22 X 10’' 

-317.52 

1.35 

14 

C0-\-0 + M^C02 + M 

3.01 X 10^^ 

1515.44 

0.0 

15 

CHO + H ^C 0 + H 2 

7.23 X 10^* 

0.0 

0.0 

16 

CHO + 0 ^ CO + OF 

3.00 X 10^2 

0.0 

0.0 

17 

CHO + OF ^ CO + F 2 O 

1.00 X 10^^ 

0.0 

0.0 

18 

CFO + O 2 ^ CO + FO 2 

4.20 X 10^2 

0.0 

0.0 

19 

CHO + M — CO + H + M n 

1.86 X 10^'^ 

8551.42 

- 1.0 

20 

CH2O + F ^ CFO + F 2 

1.26 X 10» 

1094.49 

1.62 

21 

CH2O + 0 ^ CHO + OH 

3.50 X 10^3 

1768.01 

0.0 

22 

CH2O + OF ^ CFO + F 2 O 

7.23 X 103 

-488.31 

2.46 

23 

CF 2 O + O 2 — CFO + FO 2 

1.00 X 10^^ 

20085.61 

0.0 

24 

CF 2 O + CF3 ^ CHO + CF4 

8.91 X 10-^3 

-481.09 

7.4 

25 

CH2O + M^ CHO + H + M 

5.00 X 10^® 

38487.40 

0.0 

26 

CHz + 0 ^ CH2O + F 

8.43 X 10^3 

0.0 

0.0 

27 

CF 3 + OF ^ CH2O + F 2 

8.00 X 10^3 

0.0 

0.0 

28 

CFs + O 2 - CF 3 O + 0 

4.30 X 10^3 

15503.20 

0.0 

29 

CF 3 + O 2 - CH 2 O + OF 

5.20 X 10^3 

17559.87 

0.0 

30 

CF 3 + FO 2 ^ CF 3 O + OF 

2.28 X 10^3 

0.0 

0.0 

31 

CF 3 + CHO ^ CF 4 + CO 

3.20 X lO^i 

0.0 

0.5 

32 

CHz + CHz - C2F5 + F 

4.90 X 1Q33 

5905.41 

0.0 

33 

CF 4 (+M) ^ CF 3 + F(+M) 
K) 

1.19 X 1033 

53829.68 

-4.911 


Kf 

7.05 X IQi® 

52788.95 

-0.558 


a, 6, c 

0.555 

405.62 

4580.87 

34 

CF 4 + F ^ CF 3 + F 2 

7.80 X 10® 

3896.85 

2.11 

35 

CF 4 + 0 ^ CF 3 + OF 

1.90 X 10« 

4365.91 

1.44 

36 

CF 4 + 02 ^ CF 3 + FO 2 

5.60 X 10^3 

28179.99 

0.0 
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Table 2: continued 



Reaction 

A 

0 

b 

37 

OF4 -h OF ^ OF3 -1- FaO 

1.50 X 10® 

1226.79 

2.13 

38 

OF4 FOa OF3 -H FaOa 

4.60 X 10^2 

9056.57 

0.0 

39 

OF3O + ff^ OFaO -H Fa 

2.00 X 10^3 

0.0 

0.0 

40 

OF3O -1- OF ^ OFaO -h FaO 

5.00 X 10^2 

0.0 

0.0 

41 

OF3O -1- O2 ^ OFaO -1- FOa 

4.28 X 10-^"* 

-1775.23 

7.60 

42 

OF3O + M ^ OFaO -b F -1- M 

1.00 X 10^^ 

12628.68 

0.0 

43 

O2F3 -b O2 ^ OFaO -b CEO 

3.98 X 10^2 

-120.27 

0.0 

44 

O2F4 -b F ^ O2F3 -b Fa 

3.16 X lO^i 

4029.15 

0.70 

45 

O2F4 -b OF ^ O2F3 -b FaO 

3.00 X 1013 

1503.41 

0.0 

46 

O2F5 -b F ^ O2F4 -b Fa 

3.00 X 10^3 

0.0 

0.0 

47 

OaFs -b Oa ^ O2F4 -b FOa 

2.00 X lO^'-* 

2513.71 

0.0 

48 

OaF5(-bM) ^ O2F4 -b F(+M) 
Kj 

6.24 X 10*® 

21384.56 

-6.80 


Kf 

4.97 X lOi® 

18549.48 

0.73 


a, &, c 

0.667 

653.88 

8733.74 

49 

CaF6(-bM) ^ OF3 -b OF3(-bM) 
Kj 

2.23 X 10®^ 

49895.06 

-11.992 


Kf 

7.10 X 10^® 

46866.83 

-2.792 


a, 6, c 

0.805 

302.71 

10730.56 

50 

OaFe -b F O2F5 + Fa 

5.40 X 10^ 

2621.95 

3.50 

51 

OaFe -b OF ^ O2F5 -b FaO 

2.20 X 10*^ 

565.28 

1.90 

52 

OaFe -b OF3 ^ O2F5 -b OF4 

5.50 X 10-^ 

4173.48 

4.00 


tUnits are in moles, seconds, centimeters, and Kelvins 
Third-body efficiencies 


(6) H2 = 1 . 9 , 02 = 2.6, N2 = 2.6, H2O = 9 . 5 , CO = 2.6, CO2 = 2.6 
(12) H2 = 2.9,02 = 1.2, iVa = 1 . 2 ,.ff 20 = 18 . 5 , CO = 2 . 1,002 = 4.3 
( 14 ) O2 = 12 . 0 , iVa = 2 . 0,00 = 3 . 0 , 00 a = 7.0 
( 19 ) Fa = 1 . 87 , FaO = 8.12 

( 25 ) Fa = 2 . 9 , 0 a = 1.2, Fa = 1.2, FaO = 18 . 5,00 = 2 . 1 , OO2 = 4.3 
( 42 ) Fa = 2 . 9 , 02 = 1.2, Fa = 1.2, FaO = 18 . 5,00 = 2.1, OO2 = 4.3 
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